Please note this document is a working draft.

Introduction

Background

Cities have used urban green spaces as a tool to enhance the health of residents, increase biodiversity, and mitigate the effects of climate change. Research has shown that urban green space can alleviate urban-heat-island effects, improve air quality, and encourage outdoor socialization and physical activity.1 However, urban green spaces are more frequently supported in higher-income, predominantly white neighborhoods, exacerbating inequity in health.2 Several U.S. cities, including Denver, have proposed goals to equitably increase green space in their cities.

Often, urban planners and residents use turf and exotic plants to expand green space, which are are not ideal for pollinator health and demand high amounts of water, which may not be sustainable in the mountain west, a region expected to experience more severe drought due to climate change. While a growing body of research suggests exposure to urban green space can improve human health and prevent premature mortality,1,3 little is known about the potential impact of native vegetation on human health in cities in the U.S. mountain west, a region with a growing population vulnerable to the impacts of climate change.

Objective

We aim to estimate the number of premature deaths that would be prevented by the implementation of various policy scenarios adding green space, including native vegetation, to the City of Denver. These policy scenarios were informed by conversations with local stakeholders in advocacy, research, and governance.

Methods

Measurement of existing green space and definition of native-plants

We measured existing green space in the City of Denver by the normalized difference vegetation index (NDVI) measured by the Landsat-8 satellite at a spatial resolution of 30 square meters. To define native-plants greenness, we measured the NDVI of several places in Denver (places listed here)comprised of various levels of native and adapted plants during cloud-free summer (May, June, July, August) days over five years (2016-2021). Over these days in these areas, the median NDVI was as low as 0.3228098 at the Denver Botanic Gardens Green Roof, which is a mix of native and other plant types and exists in an urban environment and as high as 0.704 in City Park, which is 0% native and consists of more watered turf. A plot in Denver Botanic Gardens that is 100% native plants had a median NDVI of 0.49.

To acknowledge the lack of precision in these NDVI estimates and the range of NDVI levels that may correspond to a native-plants intervention, we defined NDVI using 0.49 (Denver Botanic Gardens 100% native plot) as the primary analysis and 0.322 (Denver Botanic Gardens Green Roof) as an alternative value.

Stakeholder engagement to develop policy scenarios

We conducted ten interviews with local stakeholders who are actively engaged in advocacy, policy, and governance related to native plants and green space in Greater Denver. Following these interviews, we developed four policy scenarios informed by the collective set of visions and priorities articulated by the stakeholders. The scenarios range from short-term goals to more ambitious, perhaps less immediately realistic visions. The scenarios fall into five categories.

  1. Homogenous native-plants greening: Under the first set of scenarios, we consider adding native plants homogenously across all block groups, without respect to the specific location. Specifically, we consider two scenarios:

    • Set 30% of the block group’s area to be as green as native plants, corresponding with a goal set by some scientists and advocatess to protect 30% of lands and oceans by 2030[@dinerstein2019];

    • Set 20% of the block group’s area to be as green as native plants, a level that, according to some of our interviewees is sufficient to support pollinator health.

  2. Equity-driven native plants greening: Under the second set of scenarios, we propose adding native plants to the same amount of area as above (30% or 20% of Denver’s area) but rather than applying the intervention homogeneously across all black groups, we prioritize application of the native plants intervention to those areas of greatest need from a health-equity perspective. Both the State of Colorado and the City of Denver have created equity indices to guide their work. The Colorado Department of Public Health & Environment (CDPHE) describes their equity index here4 and makes the index available for download here.5 The City of Denver Department of Public Health & Environment (DDPHE) describes the neighborhood-level equity index here6 and makes the data available for download here.7 We use both definitions to consider a total of four scenarios:

    • Set 30% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the CDPHE definition of health equity.

    • Set 20% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the CDPHE definition of health equity.

    • Set 30% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the DDPHE definition of health equity.

    • Set 20% of Denver’s area to be as green as native plants, and apply this intervention to those block groups most in need based on the DDPHE definition of health equity.

  3. Create native buffers around riparian areas (rivers, streams, lakes, and ponds) of the following sizes:

    • 200 feet (ideal for ecosystem health but possibly not realistic)

    • 100 feet (compromise)

    • 50 feet (most realistic; less good for ecosystem health)

  4. Initiatives related to green infrastructure and stormwater management. (Please refer to corresponding section for details.)

  5. Replace a portion of parking-lot surface with native plants:

    • 100% of the existing parking area (unrealistic, out of curiosity)
    • 50% of the existing parking area (perhaps realistic)
    • 20% of the existing parking area (most realistic)

Visualization of scenarios

This section elaborates on the scenarios and visualizes them with maps.

Scenario 1: add native plants to all census block groups

Under the first scenario, which we view as the most ambitious, we consider the expected impact on mortality by applying a green-space intervention to all census block groups, without specifying where. This scenario was informed by conversations with colleagues who stated that

We first measured the mean NDVI of each census block group on July 4, 2021. The weighted block-group-level mean is presented here, weighted by the proportion each 30 square-meter pixel covers by the census block group. For example, if half of a pixel overlaps the block group, it receives a weight of 0.5 in the weighted average. We also removed bodies of water before measuring NDVI.

For each definition of native-plants NDVI, We excluded census blocks whose baseline NDVI was above the threshold. This map visualizes the NDVI threshold at 0.4.

And this one defines the native-plants threshold at 0.5.

Scenario 2: add native plants prioritizing census block groups with most health-equity need

Under this set of scenarios, we assumed we had a “budget” of greening of either 30% or 20% of Denver’s area, and we applied that budget to those with most need from a health-equity standpoint based on the CDPHE and DDPHE definitions, described above. As the below maps show, the two health-equity definitions from the State and from the City overlap considerably.

30% of area prioritizing block groups of greatest need according to two health-equity definitions

20% of area prioritizing block groups of greatest need according to two health-equity definitions

Scenario 3: add native-plants buffers to riparian areas

We measured NDVI in a 200-foot buffer, a 100-foot buffer, and a 50-foot buffer around all bodies of water in Denver. We downloaded bodies of water from OpenStreetMap (code here). We defined residential exposure to these riparian areas as those individuals living within a 500-meter buffer, following research on green-space and health that have defined green-space exposure based on residential proximity.3 We estimated the number of people in this buffer by multiplying the population density of the census block group by the intersecting area.

The below map depicts mean NDVI in the portions of census block groups that intersect a 200-foot buffer as well as those pieces that intersect the part of the 500 m buffer that would not be intervened upon, i.e., the part between the 200-foot buffer and the edge of the 500 m buffer.

Scenario 4: green infrastructure

We spoke with representatives at the local Office of Green Infrastructure who described three categories of initiatives–some planned, some aspirational–that could add native or adapted plants to the city.

  • Large stormwater retention projects:
    • Large ponds or basin located on public property that treat collect and treat stormwater after it has been collected in a storm pipe. They are usually vegetated, often with native or adapted plants.
    • Our colleague expected about 75% of the facility’s footprint would have native or adapted plants.
  • Green streets; according to our colleague:
    • Historically, about 2.7 miles of green streets each year, and each mile equates to about 0.15 acres of native or vegetated landscape.
    • Short term goal: increase output to 5.0 miles of green streets; same amount of vegetation per green mile
    • Aspirational goal: output to 5.0 miles; increase vegetated area to 0.75 acres per green mile
  • Proposed stormwater controls on new or re-development:
    • Our colleague stated that as properties develop or redevelop, they may be required to include stormwater runoff control measures to offset negative impacts to flooding and water quality downstream of the site associated with the impervious surfaces added during the development. These stormwater control measures are often green-on-the-ground practices vegetated with native or adapted plants.
    • The requirements may differ by parcel size. Our colleague estimated the following number of sites that may require stormwater control by parcel size:
      • greater than 1 acre: about 100 sites per year
      • 0.5-1.0 acre: about 25 sites per year
      • less than 0.5 acre: 400 sites per year

4.1. Large stormwater retention ponds

We were provided a list of planned projects throughout Denver. Per our conversations, we assume that about 75% of the project’s footprint would consist of native or vegetation. Like for riparian areas, we defined exposure to the projects as those individuals living within a 500 m buffer of the projects. A map of baseline NDVI of the projects themselves and of a 500 m buffer around the projects is below.

4.2. Green streets

Not yet complete.

4.3. Site new or re-development storwmater controls

The approach for estimating the health impact of the redevelopment controls differs from scenarios 2 and 3.1. because we will not know where the re-developments will occur. We thus simulated possible locations. Per our conversations, we anticipate the following number of parcels will be subject to these rules per year:

  • greater than 1 acre: about 100 sites per year
  • 0.5-1.0 acre: about 25 sites per year
  • less than 0.5 acre: 400 sites per year

We gathered data on parcels from Denver’s Open Data Portal (Existing Landuse 2018) and measured their area. A subset of these parcels near Union Station is mapped below.

Then, from all parcels, we sampled 100 sites of size greater than 1 acre, 25 sites of size 0.5-1 acre, and 400 sites less than 0.5 acres. One such sample appears below.

From this point, we followed the same framework as for other scenarios. The baseline NDVI of the sampled parcels and that of their corresponding 500m buffers is visualized below.

Scenario 5: add native plants to parking lots

Finally, we propose replacing a portion of parking-lot surface with native plants. Denver’s area is about 9% parking lot, so even replacing a small amount of parking-lot surface could have a large impact on the total area of urban greening. We obtained spatial data on parking lots](https://www.denvergov.org/opendata/dataset/parking-lots-2016)) from the City of Denver Open Data Catalog.

Total area (mi2), parking lots Total area (mi2), Denver proportion parking lots
13.57 155 0.09

We similarly measured the NDVI on the parking lots as well as the NDVI in the census block groups within a 500-m buffer radius of any parking lot (which is almost the whole city). We present a small subset of these areas near Union Station:

Health-impact assessment methods

We estimated the number of premature deaths averted under each scenario by following a recent meta-analysis that estimated that for every 0.1 unit increase in exposure to NDVI, the relative risk of premature death decreases about 4% (pooled risk ratio of 0.96; 95% confidence interval [CI]: 0.94, 0.97).3 We used this risk ratio to estimate the population attributable fraction8 corresponding to the proposed changes in NDVI from the baseline levels for each scenario to the alternative native-plants level.

Under the first two scenarios, the areal unit is the block group itself, and because we do not specify where, exactly, in the block group the greening intervention would occur, we define NDVI exposure as the average NDVI in the block group and compare that value with the alternative native-plants value. We use 5-year American Community Survey Data (2016-2020) to estimate the population in 5-year age groups in each census block group. For the other scenarios, we draw a 500-m buffer around the intervention area, as shown above. To estimate the population affected, we multiply the population density of the block group in that age group by the area of the block group covered by the intervention area.

To estimate the number of premature deaths prevented in each age group and area, we multiply the population-attributable fraction by the baseline mortality rate by the population size. We restrict analyses to adults aged 30 and above (subject to change) following the age range of many of the cohort studies reviewed3 and exclude northeastern census block groups near the airport (subject to change; perhaps not for all scenarios).

We estimate uncertainty in both the block-group-level population estimates and the relative risk estimate (95% CI: 0.94, 0.97) by resampling from normal distributions and taking the 2.5th and 97.5th percentiles over 1,000 replications. This code has additional detail.

Results

Results of each scenario appear in the table below. The first table describes, for each scenario, the population affected, the cumulative area of the treatment itself (e.g., the area of the retention ponds or parking lots) and the cumulative area of the 500 m residential buffer surrounding the treatment area. Note for Scenarios 1 and 2, the treatment area and the residential-buffer area are the same (the area of the affected census block groups) because we do not say exactly where in the block groups the greening treatment would occur.

Scenario Sub-scenario NDVI, native Area, treatment (mi2) Area, residential buffer (mi2) Pop. affected, est. Pop. affected, 95% LL Pop. affected, 95% UL
all-bg 100-pct 0.32 1,053 1,053 128,501 125,014 131,828
all-bg 100-pct 0.49 2,430 2,430 369,667 363,881 376,485
all-bg 30-pct 0.32 1,053 1,053 128,501 126,536 133,894
all-bg 30-pct 0.49 2,430 2,430 369,667 364,577 376,402
all-bg 20-pct 0.32 1,053 1,053 128,501 125,874 131,728
all-bg 20-pct 0.49 2,430 2,430 369,667 367,377 374,741
all-bg-equity cdphe-def-30-pct 0.32 486 486 45,980 43,849 46,838
all-bg-equity cdphe-def-30-pct 0.49 792 792 104,970 100,354 107,172
all-bg-equity cdphe-def-20-pct 0.32 371 371 29,417 27,778 30,202
all-bg-equity cdphe-def-20-pct 0.49 538 538 61,145 58,681 62,439
all-bg-equity denver-def-30-pct 0.32 454 454 39,260 35,946 40,533
all-bg-equity denver-def-30-pct 0.49 792 792 107,015 101,182 110,362
all-bg-equity denver-def-20-pct 0.32 374 374 27,246 25,597 28,212
all-bg-equity denver-def-20-pct 0.49 538 538 60,046 56,661 62,047
riparian 200-ft 0.32 68 598 60,816 58,486 62,807
riparian 200-ft 0.49 147 1,267 174,284 171,598 176,859
riparian 100-ft 0.32 34 598 60,816 59,115 63,436
riparian 100-ft 0.49 72 1,267 174,284 170,346 176,380
riparian 50-ft 0.32 17 598 60,816 59,379 63,371
riparian 50-ft 0.49 36 1,267 174,284 172,212 177,583
ogi ogi_proj 0.32 9 232 40,046 37,609 42,039
ogi ogi_proj 0.49 17 440 74,882 72,539 77,567
ogi parcel 0.32 1 786 102,819 100,632 107,018
ogi parcel 0.49 2 1,911 295,660 290,277 300,392
prkng 100-pct 0.32 17 1,053 129,060 125,226 132,507
prkng 100-pct 0.49 24 2,423 363,994 355,196 369,639
prkng 50-pct 0.32 17 1,053 129,060 125,537 131,053
prkng 50-pct 0.49 24 2,423 363,994 361,546 370,440
prkng 20-pct 0.32 17 1,053 129,060 125,119 131,572
prkng 20-pct 0.49 24 2,423 363,994 359,064 366,026

The second table shows, for each scenario, the estimated number of deaths prevented and the number of deaths prevented per 100,000 people affected.

Scenario Sub-scenario NDVI, native Deaths prevented, est. Deaths prevented, 95% LL Deaths prevented, 95% UL Deaths prevented per 100k, est. Deaths prevented per 100k, 95% LL Deaths prevented per 100k, 95% UL
all-bg 100-pct 0.32 40 36 46 31 28 35
all-bg 100-pct 0.49 219 206 239 59 56 65
all-bg 30-pct 0.32 12 10 15 9 8 11
all-bg 30-pct 0.49 64 59 70 17 16 19
all-bg 20-pct 0.32 8 7 9 6 5 7
all-bg 20-pct 0.49 42 41 46 11 11 12
all-bg-equity cdphe-def-30-pct 0.32 4 3 5 9 8 10
all-bg-equity cdphe-def-30-pct 0.49 22 21 24 21 20 22
all-bg-equity cdphe-def-20-pct 0.32 2 1 2 6 5 7
all-bg-equity cdphe-def-20-pct 0.49 8 8 10 14 13 15
all-bg-equity denver-def-30-pct 0.32 3 2 3 7 6 8
all-bg-equity denver-def-30-pct 0.49 19 17 19 17 16 18
all-bg-equity denver-def-20-pct 0.32 1 1 2 5 4 6
all-bg-equity denver-def-20-pct 0.49 7 6 8 12 11 13
riparian 200-ft 0.32 1 0 1 1 1 2
riparian 200-ft 0.49 7 6 8 4 3 5
riparian 100-ft 0.32 0 0 0 0 0 0
riparian 100-ft 0.49 2 2 2 1 1 1
riparian 50-ft 0.32 0 0 0 0 0 0
riparian 50-ft 0.49 1 0 1 0 0 1
ogi ogi_proj 0.32 0 0 1 1 0 2
ogi ogi_proj 0.49 1 0 2 2 1 3
ogi parcel 0.32 0 0 1 0 0 1
ogi parcel 0.49 1 0 1 0 0 0
prkng 100-pct 0.32 12 10 13 9 8 10
prkng 100-pct 0.49 40 37 43 11 10 12
prkng 50-pct 0.32 6 5 7 5 4 5
prkng 50-pct 0.49 20 19 21 6 5 6
prkng 20-pct 0.32 2 2 3 2 2 2
prkng 20-pct 0.49 8 7 8 2 2 2

Next steps

As of April 27, 2022:

  • Add green streets HIA

  • Add economic-impact attributable to health

References

1.
2.
Wen M, Zhang X, Harris CD, Holt JB, Croft JB. Spatial Disparities in the Distribution of Parks and Green Spaces in the USA. Annals of Behavioral Medicine. 2013;45(S1):18-27. doi:10.1007/s12160-012-9426-x
3.
Rojas-Rueda D, Nieuwenhuijsen MJ, Gascon M, Perez-Leon D, Mudu P. Green spaces and mortality: A systematic review and meta-analysis of cohort studies. The Lancet Planetary Health. 2019;3(11):e469-e477. doi:10.1016/S2542-5196(19)30215-3
4.
CDPHE. Climate Equity Data Viewer. Published online September 13, 2021. https://storymaps.arcgis.com/stories/46bf289f92bc4629a0a1266de4bb7f97
5.
6.
7.
Denver. Denver open data catalog: Equity index 2020 - neighborhood. https://www.denvergov.org/opendata/dataset/city-and-county-of-denver-equity-index-2020-neighborhood
8.
Murray CJ, Ezzati M, Lopez AD, Rodgers A, Vander Hoorn S. Comparative quantification of health risks: Conceptual framework and methodological issues. Population Health Metrics. 2003;1(1):1. doi:10.1186/1478-7954-1-1


Copyright © 2022 Michael D. Garber